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ABSTRACT 

We present the discovery of [OIII] 5007 emission associated with the black hole X-ray binary recently 
identified in a globular cluster in the Virgo elliptical galaxy NGC 4472. This object is the first con- 
firmed black-hole X-ray binary in a globular cluster. The identification of [OIII] 5007 emission from the 
' black-hole hosting globular cluster is based on two independent fiber spectra obtained at the VLT with 

I FLAMES, which cover a wavelength range of 5000-5800 A at a spectral resolution of about 6000. In 

_ ■ each of these spectra we find an emission line at 5031.2 A with an uncertainty of several tenths of an 

angstrom. These are consistent with [OIII] 5007 emission at the 1475 ± 7 kms~^ radial velocity of the 
, globular cluster previously determined from an analysis of its absorption lines. This agreement within 

the small uncertainties argues strongly in favor of the interpretation of the line as [OIII] 5007 emission 
from the black-hole hosting globular cluster. We also find that the emission line most likely has a velocity 
" ' . width of several hundred kms^^. Such a velocity width rules out a planetary nebula explanation for 

the [OIII] 5007 emission and implicates the black hole as the source of the power driving the nebular 
I emission. 

Oh' Subject headings: galaxies individual (NGC 4472) — galaxies: star clusters — globular clusters: 

Q , general — X-rays: binaries — X-rays: galaxies: clusters 

H ■ 

■ 1. INTRODUCTION black-hole X-ray binary in a globular cluster. As such it 

Maccarone et al. (2007, hereafter MKZR07) presented °P^,^' .^P new possibilities to address questions about the 

,1 J- r ui 1 1 1 V u- ■ 1 u 1 evolution of black holes m dense stellar systems. It also 

. the discovery oi a black-hole X-ray binary m a globular . , . i 

^ , cluster in the Virgo elliptical NGC 4472. This is the first f ' ^ unique opportunity to study m a nearby system 

>, • ui 1 u 1 V u- • 1 u 1 1 how the energy from a highly luminous accreting black 

unambiguous black- hole X-ray binary m any globular cms- i i • n i i i • , •, t t • i 

' 4. rrvT- li • u J inn 1 vn^nj at j- t, holc IS icd back into its surrounding medium m a dense 

m , ter. This result is based on a 100 ks XMM-Newton ob- , , mi- • i • ; , , • , i • 

O.. r .1 -ir- 11- 1 AT/-</--i AATr, J IT. stcllar systcm. io begin considering these questions, this 

■ servation of the Virgo elliptical NGC 4472 and its glob- , , • •,• , ;• , , r ;i • , , , 

^ ' ular clu^ter^ STiecificflllv MKZR07 identified an X rav Presents initial optical spectroscopy of this black- 

^ uiai ciUoLeiD. o ueciiicaiiy, iviivzjixu/ loeiiLiiieci an -A.-iay ^ ^ -^r ^ • i ii i i i i i • i ■ i •, • 

rf^ ' . , . c J 1 1. 1 1 4- r hole X-ray binary and the globular cluster m which it is 

; . source m a spectroscopically confirmed globular cluster ot ^ , , , , . . „ , . „ , ., 

, -rr. 1,- ,. , 1 T\jr<r' AA^r, -li V 1 • found. We present the observations m Section 2, describe 

0> ' the Virgo elhptical galaxy NGC 4472 with an X-ray lumi- ,, r .^ ■ ^ ■ ■ o i i- ;i • 

O., r . in3q -i-c i-Ui-u u J i- the results oi their analysis m section 3, and discuss their 

nosity 01 4 X 10 ergss lor which the observed count . ■ , 

X — J. -JT, fj- r c \. rpii-i implications m Section 4. 

l ■ rate varied by a factor oi seven over a tew hours. Ihe high ^ 

X-ray luminosity, about ten times I/Edd for a neutron star, 

rules out any single object other than a black hole in the 2- observations 

old steUar population of the globular cluster (MKZR07). As part of a study of the kinematics of the globular 

, However, a high X-ray luminosity alone does not guaran- cluster system of Virgo elliptical galaxy NGC 4472, we 

. tee a black hole source in extragalactic globular clusters, used the FLAMES multi- fiber spectrograph (Pasquini et 

as multiple neutron stars are a possibility in these sys- al. 2004) on the VLT/UT2 to obtain spectra of globular 

tems which are unresolved in X-ray data. Variability has clusters over a wide area around NGC 4472 (Zepf et al. 

long been recognized as a key for distinguishing a black 2007). The globular cluster targets were selected from 

hole X-ray binary from multiple neutron stars in extra- the multicolor, wide-field imaging study of Rhode & Zepf 

galactic globular clusters (e.g. Kalogera, King, & Rasio (2001, hereafter RZOl). RZOl used BVR MOSIAC im- 

2004). Therefore, the clear variability observed during the ages from the Mayall 4-m to identify candidate globular 

XMM observation unambiguously indicates that the high clusters over an approximately 35' x 35' field of view. The 

X-ray luminosity in this source originates from a black hole spectroscopic data were obtained with FLAMES in its GL 

(MKZR07). RAFFE/MEDUSA mode, which offers 130 fibers over a 25' 

This object is the first and currently only confirmed diameter field of view. By observing through four fiber se- 
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tups, each slightly offset in a different direction from the 
center of the galaxy, we were able to obtain spectra for 
nearly 400 globular clusters over a large range of distances 
from the center of the galaxy, with good azimuthal cov- 
erage, and many overlaps between pointings to provide 
checks on the accuracy of the radial velocities of the clus- 
ters. 

The observations were carried out in service mode at the 
VLT/UT2 with two fiber setups observed in Jan 2004 and 
two observed in Apr 2005. All of the data were obtained 
with grating LR4, which covers the range 5015 A to 5831 A 
with a resolving power A/A(A) of 6000. This range covers 
many stellar absorption features including strong Mg and 
Fe lines and is therefore excellent for determining the ra- 
dial velocity of globular clusters. The fairly high spectral 
resolution also results in small velocity uncertainties (see 
Bergond et al. 2006, Zepf et al. 2007). For studying pos- 
sible emission lines from the black-hole X-ray binary, the 
narrow wavelength range is imperfect. However, it does 
provide a first look at the spectrum of the unique source, 
and the spectral resolution corresponding to a velocity 
width of about 50 kms~^ is useful for assessing whether 
any of the lines are broader than this value. 

The bright X-ray source XMMU 122939.7-^075333 iden- 
tified in the XMM and Chandra images is located at 
the position of the globular cluster candidate RZ 2109 in 
the Rhode & Zepf 2001 cluster sample. The astromet- 
ric matching of these catalogs was described in MKZR07, 
with the agreement good to < 0.4", a level similar to the 
positional uncertainties in the catalogs themselves. As 
discussed in MKZR07, this excellent agreement gives a 
negligible probability (< 10^®) that the X-ray source and 
globular cluster are not associated with one another. This 
globular cluster was observed in two of the FLAMES fiber 
setups described above. One of these was observed for a 
total of 4.3 hours in Jan 04, and the other for 3.2 hours in 
Apr 05. The reduction of the FLAMES multifiber spectra 
was carried out in the same way as our previously pub- 
lished analysis of similar data for globular clusters around 
NGC 3379 (Bergond et al. 2006). 

Our standard analysis of the radial velocity of RZ 2109 
from the cross-correlation of its absorption lines with stel- 
lar templates gives a radial velocity of 1477 ± 7 kms~^ for 
the Jan 04 observations, and 1460 ± 19 kms~^ for the Apr 
05 observations. We adopt a weighted average of 1475 ± 7 
kms^^as the best estimate of the globular cluster radial 
velocity (in MKZR07 we used the 1477 kms"^ of the ob- 
servation with the smallest error). The agreement between 
the two gives high confidence in this radial velocity, and 
its value unambiguously identifies the object as member of 
the NGC 4472 globular cluster system, which has a mean 
radial velocity of 1018 kms~^ and a velocity dispersion of 
about 300 kms-i (e.g. Zepf et al. 2000). 

3. EVIDENCE FOR [OIIl] 5007 EMISSION 

To test for the presence of emission lines in the opti- 
cal spectrum of XMMU 122939.7-^075333, we visually in- 
spected the two independent FLAMES spectra of the ob- 
ject described above. As shown in Figure 1, both spectra 
show evidence for an emission line at a wavelength corre- 
sponding to [OIII] 5007 at the radial velocity observed for 
the host globular cluster. Although the line is only seen 



at modest signal-to-noise in each spectrum, the fact that 
it is observed in the same place in the two independent 
spectra, and that it is located precisely at the wavelength 
of [OIII] 5007 redshifted by the radial velocity of the clus- 
ter strongly supports this identification of the line. More 
specifically, we find the line center to be 5031.5 A in Jan 
04 spectrum, while in Apr 05 spectrum, we find a line 
center of 5030.8 A with an uncertainty of several-tenths 
of an angstrom, based on both gaussian fits and simple 
mean wavelengths of the continuum subtracted emission 
line. We adopt 5031.2 ± 0.3 A as the best estimate of the 
emission line center. Based on a zero redshift wavelength 
of 5006.8 A for the [OIII] line, we then find a radial veloc- 
ity of the line of about 1462 kms^^ with an uncertainty of 
about 18 kms^^, in excellent agrecnicnt with 1475 kms~^ 
found for the stellar absorption lines of the host globular 
cluster. 

A second key feature of the line seen in each spectrum 
is that it appears to be broader than the spectral reso- 
lution of the data. Specifically, standard gaussian fits to 
the line give a FWHM of 4 — 6 A in both independent 
spectra. This is much larger than the 0.9 A (4.5 pixels) 
FWHM of the effective instrumental resolution. The [OIII] 
5007 line then is clearly intrinsically broad, with a width 
of about 200 — 350 kms~^. This measured width does not 
rule out either a weaker, broader component, or detailed 
structure within the line. However, it does clearly indicate 
the overall [OIII] 5007 emission is spread over at least a 
few hundred kms~^. 

This observed velocity width strongly implicates the 
black hole as the source of the energy driving the [OIII] 
5007 emission. Without a constraint on the line width, 
the possibility would remain that a planetary nebula in 
the globular cluster is the source of the [OIII] emission. 
While these are rare, there are two known cases among the 
hundreds of extragalactic globular clusters with published 
spectroscopy (Bergond et al. 2006 and Pierce et al. 2006 
for a NGC 3379 globular cluster, and Minniti & Rejkuba 
2002 for a NGC 5128 globular cluster). However, plane- 
tary nebulae have velocity widths of tens of kms~^ (Acker 
et al. 1992 and references therein), and therefore the ob- 
served velocity width of several hundred km eliminates 
the possibility that the optical emission is from a PNe. 

The luminosity of the emission line can also be esti- 
mated from the spectra. Although flux standards were 
not observed as part of this program, the target globular 
cluster has a known magnitude of V=21.0 (RZOl). As a 
result, we can fix the flux in the continuum to be that 
for a globular cluster with V=21.0, and use that flux cal- 
ibration to estimate the observed flux in the [OIII] 5007 
line. We can also check that the detailed flux distribution 
across the V band agrees with that expected from mod- 
els of stellar populations (e.g. Bruzual & Chariot 2003) 
like that of the globular cluster and its known colors. It 
is worth noting that the uncertainty in the determination 
of the total counts in the omission line is substantial, as 
described below, and thus the uncertainty in this fltix cal- 
ibration contributes negligibly to the error budget of the 
line luminosity determination. 

It is in this way that the flux calibration was estimated 
for the two spectra shown in Figure 1. An independent 
check on this calibration can be made by comparing the 
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observed signal to noise to that expected given the cluster's 
known magnitude and the observing conditions, based on 
the well characterized throughput of the FLAMES instru- 
ment. We find good agreement between these, providing 
further evidence that the flux estimate is reasonable, and 
that the uncertainty in the line luminosity is dominated 
by the measurement of the line itself. 

The [OIII] line luminosity can then be determined by 
measuring the total line flux in the calibrated spectrum, 
and adopting a distance to NGC 4472 of 16 Mpc (Macri et 
al. 1999). The most uncertain part of this procedure is de- 
termining the total fux in the line, and in particular setting 
the level of the continuum and the range of wavelengths 
over which to determine the line flux. Our best estimate 
from both of our spectra in Figure 1 is that L([OIII] 5007) 
is about 4 x 10^^ ergss"^ . However, if the line is broader 
than we are able to clearly discern in our data, we esti- 
mate this value could be roughly twice as large and we 
also find it might be possible to to fit the line with a flux 
that is smaller by 50%. Therefore, we attribute a factor of 
two uncertainty to our L([OIII] 5007) ~ 4 x 10^^ ergss"^ 
estimate. 

4. DISCUSSION 

We have shown that optical spectroscopy reveals [OIII] 
5007 emission associated with the black-hole X-ray binary 
system in a globular cluster around the elliptical galaxy 
NGC 4472. This [OIII] 5007 emission clearly estabhshes an 
interaction between the X-ray binary and material around 
the black hole system. As such it also provides a poten- 
tially valuable testbed for understanding the relation be- 
tween the binary evolution and accretion history of black 
holes and their interaction with the surrounding interstel- 
lar medium. Our spectra have also been able to establish 
the general velocity width of the [OIII] 5007 line. The sev- 
eral hundred kms~^ width of the line both demonstrates 
the emission is associated with the black-hole X-ray bi- 
nary since no other cluster source could produce it, and 
also shows the material may be driven out of the globu- 
lar, as the velocity width is much larger than the escape 
velocity of globular clusters. 

[OIII] 5007 can be produced in principle by either col- 
lisional or photoionization. We consider each of these 
possibilities and their implications for the mass of the 
black hole. In the collisional case, a strong wind is gen- 
erated by the black-hole X-ray binary, which then drives 
a shock front into the interstellar medium of the globular 
cluster. We can write standard solutions for the radius 
R of this shock front at time t as R{t) = 28[(M/10-*^ 
M0/yr)(i;o/2OOOO kms-^)^{0.1cm-^ /no)]^/^{t/10^yv)^/^ 
pc (e.g. Castor, McCray, & Weaver 1975), and the cur- 
rent velocity as V{t) = 3/5R{t)/t, where M is the mass 
loss rate in the wind, vq the initial velocity, and no the 
density of the interstellar medium. We now consider what 
likely values for no, M, vq, and t are and how these might 
be constrained by current and future observations. 

The density of the interstellar medium of globular clus- 
ters, no, is only well-known in a few cases. The best con- 
straint is probably for 47 Tuc, for which the gradient in 
pulsar dispersion measurements with radial position in the 
cluster gives a free-electron density of 0.067 ± 0.015 cm~^ 
(Preire et al. 2003). Similar studies of pulsar dispersion 



measurements give a similar or slightly higher density in 
M15 (Freire et al. 2001), and a recent study may have de- 
tected somewhat more neutral hydrogen in M15 (van Loon 
et al. 2006). Most other studies of globular clusters have 
only been able to place upper limits which are typically 
well above the adopted value. That there would be some 
interstellar medium in a globular cluster seems natural, as 
mass loss from the cluster stars provides a continual source 
of interstellar material. Even if many processes cause clus- 
ters to also losc^ ISM. the velocitic^s with which this mass 
is driven out would have to be extraordinarily high or else 
the steady state density of the ISM will be generally simi- 
lar to the fiducial value adopted here (Pfahl & Rappaport 
2001). 

The radius and current velocity of the shocked region 
also depend on the kinetic luminosity in the wind, Mvq. 
In the context of super-Eddington accretion onto a stellar 
mass black hole, this should be close to the Eddington lu- 
minosity of the source, and should be independent of the 
actual accretion rate, providing it is well above the ac- 
cretion rate require to yield the Eddington luminosity for 
an accretion efficiency of O.lc^ (Begclman, King & Pringle 
2006). The observed X-ray luminosity of = 4 x lO"'^'' 
ergs/sec and and inner disk radius of ~ 4000 km inferred 
from the X-ray data (MKZR07) can be well fit by this 
model with a 10 Mq black hole accreting at about 40 times 
the rate required to give the Eddington luminosity. 

In this case, the predicted initial wind velocity vq is 
20000 km/sec and M is 8 x lO"'' Mq. These then give a 
timescale t of about 2 x 10^ years for the bubble to ex- 
pand to the point where the one-dimensional velocity of 
the outflow is about 150 km/sec. The corresponding size 
of the bubble is about 40 pc, which is within the tidal 
radius expected for a fairly massive outer globular cluster 
like RZ 2109. Taken at face value, a model in which M 
is constant over the above timescale loses somewhat more 
mass than the 1 Mq upper limit for a donor star in the 
globular cluster. However, the detailed application of this 
model involves many uncertainties including the level of 
isotropy of the accretion flow, uncertainties in the radii in- 
ferred from X-ray spectroscopy, and inhomogeneity in the 
ISM along with the effect of input of additional material 
from ongoing stellar mass loss in the globular cluster. It 
would then seem a wind-driven system with a constant M 
remains feasible. 

Alternatively, if the system experiences periods of qui- 
escence, the total mass lost and kinetic power are less than 
calculated assuming these are constant at the above val- 
ues. This is because the timescale is the time since the 
object first became an energetic binary, while M and thus 
the kinetic luminosity are the time averaged values since 
this time, and are thus reduced directly by the duty cy- 
cle. In either case, the calculation suggests that a sce- 
nario in which a super-Eddington stellar mass black hole 
driving a strong wind through the interstellar medium 
of the globular cluster may be able to produce the ob- 
served [OIII] 5007 emission line. The X-ray data also allow 
the possibility that the accretor is an intermediate-mass 
black hole (IMBH) emitting at sub-Eddington luminosi- 
ties (MKZR07). However, a strong wind is not expected 
in the IMBH sub-Eddington case (Proga 2007 references 
therein), so if the [OIII] 5007 emission is from shocks, then 
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it would seem the black hole has a stellar mass. 

We now consider the possibility that the [OIII] 5007 
emission line is photoionized. Because the velocity width 
of the line is significantly greater than the virial velocity 
of the globular cluster, if the line is photoionized, it must 
be produced near the black hole. For a 10^ Mq black hole, 
the velocity width indicates the line must come from about 

7 X 10^^ cm away from the black hole. This may be pos- 
sible, although it may prove difficult to have enough mass 
at the right distance from the black hole to produce the 
[OIII] 5007 line luminosity at the observed velocities. For 
completeness we note that for a stellar mass black hole the 
observed [OIII] 5007 seems unlikely to be due to photoion- 
ization. The broad lines observed in Galactic black-hole 
X-ray binaries are typically permitted lines such as the 
Balmer lines and He II 4686, as might be expected due 
to the high densities at the small distances from the cen- 
tral source required to match the velocities. In contrast, 
as described above a shock is expected if the accretor is a 
stellar mass black hole and can produce roughly the ob- 
served [OIII] 5007 emission. A final alternative is it may 
be possible for a stellar mass black hole to drive a wind 
as described earlier, but for photoionization of the wind 
material by the central source to play a significant role in 
producing [OIII] 5007. 

This analysis shows that undertanding the physical pro- 
cess which produces the [OIII] 5007 emission can place 
constraints on the mass of the black hole. Determining 
the mass of the black hole has important implications be- 
yond understanding this individual system. Specifically, it 
would seem very unlikely to find a stellar mass black-hole 
X-ray binary in a globular cluster with an intermediate 
mass black hole. This is because most black holes are be- 
lieved to be either dynamically ejected (Kulkarni, Hut, & 
McMillan 1993, Sigurdsson & Hernquist 1993) or merge to 
form a IMBH (Miller & Hamilton 2002). Any surviving 
stellar mass black hole in a cluster with an IMBH would 
have to reside in the outer regions of the cluster to avoid 
either of these fates. However, a location in the outer re- 
gion of a globular cluster with its low stellar density is 
not conducive to the stellar dynamical interactions that 
make close, accreting binaries with black hole primaries in 
globular clusters. Therefore, it appears that if it can be 
shown the black hole must have a stellar mass, then this 
globular cluster would be very unlikely to host an IMBH. 
This would then demonstrate that the Mbh ~ c relation 
for larger mass stellar systems does not hold for all stellar 



systems at globular cluster masses. 

There are several future observations that will help con- 
strain the physical parameters of this interesting system. 
One of these is expanding the wavelength coverage beyond 
the narrow region currently available. This will provide 
hne ratio diagnostics such as [SII]/Ha and [0III]/H/3 to 
determine the relative importance of coUisional and pho- 
toionization processes. As described above, this will have 
immediate implications for the mass of the accreting black 
hole, as the currently observed [OIII] 5007 would seem to 
require photoionization for an IMBH, while shock ioniza- 
tion works well in the case of a stellar mass black hole. 
Wider wavelength coverage would also lead to constraints 
on the density structure in the cluster interstellar medium, 
and possibly assessing the total power and lifetime of the 
system. 

A key parameter that may be constrained by future 
observations is the radius at which the [OIII] 5007 emis- 
sion originates. We note that at the 16 Mpc distance of 
NGC 4472, 1" corresponds to 77 pc, so images with much 
higher spatial resolution than this, such as from HST, may 
be able to resolve the nebulae and determine its size if 
it originates in shocks. Correspondingly, a tight upper 
limit would suggest photoionization, as the spatial scale in 
this case is orders of magnitude smaller and unresolvable. 
Moreover, if the line emission originates in shocks, resolv- 
ing it would allow the size and velocity width to be com- 
bined to directly give the timescale, and the total kinetic 
luminosity would then follow. This would also provide a 
constraint on whether the accretion has been constant or 
intermittent, since it is the time averaged kinetic luminos- 
ity that matters for the shock radius while the timescale 
that matters is the total lifetime of the system. 
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Fig. 1. — Plots of the optical spectra of the globular cluster shown to host a black hole X-ray binary by MKZR07. Both of the two 
independent fiber spectra from FLAMES on the VLT show emission with a central wavelength of 5031.2 A with a uncertainty of several 
tenths of an angstrom, and a line width of about 4 — 6 A . The observed wavelength strongly supports the identification of this line as [OIII] 
5007 from the globular cluster, which has a radial velocity of the 1475 kms~^ determined from its stellar absorption lines. The dotted line 
denotes the wavelength expected at [OIII] 5007 at this redshift. The 4 — 6 A width of the [OIII] 5007 line is clearly broader than the 0.9 A 
spectral resolution of the data. This width corresponds to several hundred kiiis^^, much larger than any planetary nebulae, and implicates 
the black-hole X-ray binary as the source of the energy driving this emission. The plotted data have been smoothed by the instrumental 
resolution, and the calibration of the flux given on the y-axis is described in the text. 



